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Abstract 

It is shown that a quantum Levy process in a box leads to a problem involving 
topological constraints in space, and its treatment in the framework of the path 
integral formalism with the Levy measure is suggested. The eigenvalue problem 
for the infinite potential well is properly defined and solved. An analytical ex¬ 
pression for the evolution operator is obtained in the path integral presentation, 
and the path integral takes the correct limit of the local quantum mechanics 
with topological constraints. An example of the Levy process in oscillating walls 
is also considered in the adiabatic approximation. 
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1. Introduction 


The introduction of a fractional concept in quantum mechanics with motivat¬ 
ing new implementations of non-local physics leads to many technical questions 


and often needs special care. A typical example is the 


of a particle in an infinite potential well, suggested in ll 


quantum Levy flights” 


. This “sim 


pie has evoked an active discussion in the literature on how a 

non-local operator, defined on the infinite scale, acts in a finite-size range, such 
as a quantum box. The aim of the present research, related to this question, 
is an exploration of the Levy flights in boundary value problems in a finite-size 
area. Among many possible applications of the problem, of special interest is 
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the first-passage analysis [8j, which is important for the investigation of the 
transport properties of Levy glasses Q . 

In quantum mechanics, the fractional concept has been introduced by means 
of the Feynman propa gato r for non-relativistic quantum mechanics as for Brow¬ 
nian path inte gral s [l|, [lOj . Equivalence between the Wiener and the Feynman 


path integrals 


HI, established by Kac 112|, is natural, since both are Markov 


processes. As shown in the seminal papers [1|, Il0|, the appearance of the space 
fractional derivatives in the Schrodinger equation is natural and relates to the 
path integrals approach. The path integral approach for Levy stable processes, 
considered for the fractional diffusion equation, in particular, for the fractional 
Fokker-Planck equation (FFPE), has been extended to a quantum Feynman- 
Levy measure that leads to the space fractional Schrodinger equation (FSE) 

QEoi. 


The introduction of the Levy measure in quantum mechanics is based on 
the generalization of the self-consistency condition, known as the Bachelor- 
Smoluchowski-Kolmogorov chain equation (or the Einstein-Smoluchowski-Kol- 
mogorov-Chapman equation, see e.g., [l^), established for the Wiener process 
for the conditional probability W{x,t\xo,to) 


W(x,t\xf),tQ) = / W{x,t\x',t')W{x',t'\xQ,t()) dx . 


( 1 ) 


In the case of the translational symmetry, it reads IF(x,t|a;o,to) = W{x — 
Xo,t\to). Straightforward generalization of this expression by the Levy process 
is expressed through the Fourier transform 

/ OO 

g^p(a,-xo)g-iC„t|pr (2) 

-OO 

where 0 < a < 2 and Ka is a generalized diffusion coefficient Q. Using Eqs. 
m and one defines the integration of function F[a;(T)] over the generalized 
measure 


dxi... 


dxN 


dpi... 


dpN 


xF{xi, . . . , a;jy)e-d(a:^l-a:o)Pl-l-----l-(a;iV + l-a:jv)PN]g-^!'a[|pirAt-|-...-|-|pjvrAt] ^ 
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Here, At = t/N. The Levy distribution is defined by the Fox function [l| and 
in the r.h.s. of Eq. ([3]) it is presented by means of the Fourier integral with a 
stretched exponential kernel. When F[a;(r)] = exp{— F[x(r)](ir}, with V(x) 


being a potential, expression ([3]) is a generalized Feynman-Kac formula 11 m 


One should recognize that the general form of the path integral in Eg. (131) 
does not correctly describe systems with a so-called topological constraint |16j . 
In this case integrations over coordinates xi,X 2 , ■ ■ ■, xn have finite limits. As an 
example of such system, we consider fractional quantum mechanics in an infinite 
well potential. We show that the condition of the restriction of the integration 
in both the formulation of the problem in Eq. ([3]) and correspondingly, in the 
fractional Riesz derivative must be taken into account. Otherwise, the FSE for 
the infinite well potential cannot be obtained. 

2. Levy quantum mechanics in potential well 

In complete analogy with the FFPE, fractional quantum mechanics can be 
constructed from the Feynman-Kac formula and V(x) is the potential. Following 
[l| , we determine a wave function at the moment t + At by means of Eq. m 


'ip{x, t + At) = 


dp _ 
2TTh 


zDq, Ai|p|^ 


dye 




e . (4) 


In the limit At —>■ 0, one obtains 
r) /*oo POO 

ihdt'ilj{x,t) = J dy J \p\‘^e~^P^'^~y^'ilj{y,t)dp + V{x)ilj{x,t). (5) 

If one introduces the fractional Laplacian (—A)“/^ by means of the Riesz frac¬ 


tional derivative 


li 




1 


'\p\°'^{p)dp, 


( 6 ) 


where (t>{p) = is the Fourier image, the fractional Schrodinger 

equation (FSE) is obtained from Eq. ([5]) 


ihdtip = Da{—ih)°'d)^ip + V{x)'ip . 


(7) 
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One sees that for a = 2, FSE © reduces to a conventional Hamiltonian me¬ 
chanics with I ?2 = being a half inverse mass of a particle. 

It should be noted that this equation is valid only for the smooth potentials 
V{x). In particular, for the infinite potential well 


V{x) 


0 if |a;| < L 

oo if |a;| > L. 


( 8 ) 


one cannot obtain the Schrodinger equation © from Eqs. (|3]) and ([5]). A simple 
explanation of this problem relates to the absence of the expansion over V At in 
Eq. ([4]), since exp [—^AtV{y)] is an infinitely fast oscillating function for any 
arbitrary small but finite At. Note that one first performs this expansion and 
then takes the limit At = 0. To overcome this obstacle, the cutoff of the limits 
of integration over y is performed. Moreover, this is a correct formulation of the 
problem with the topological constraints, as the infinite potential well is. This 
notion is also relevant for conventional local quantum mechanics, see e.g., [l^ . 
Therefore, FSEs ([51) and ([7]) are reduced to the well defined problem of a free 
particle in a finite-size range with the FSE 

jj pL POO 

ihdt'ip = -^J dy J \p\°‘e~^P^^~'>^'>ip(y,t)dp, (9) 

which is furnished with the boundary conditions 'ijj{x = ±L) = 0. Obviously, 
the Riesz fractional derivative © can no longer be used for the potential well, 
and in addition the Fourier image 4>{p) is not appropriately defined. 

Let us perform the Fourier inversion in Eq. m over k = p/h. This yields 

27^ 7-00 27r 7_^ 


One obtains this expression with the double differentiation for 1 < a < 2, while 
for 0 < a < 1 one differentiates only once. Therefore, we have the integration 

27r|2|^-i 


\k\ 




dk = 2 \k\ '' cos kz = 


( 11 ) 


2T{v) cos(i/7r/2) ’ 

where i/ = 2 — aforl<a<2 and z/ = 1 — a for 0 < a < 1, and F(z/) = {v — 1)! 
is a gamma function. To be specific, we consider 1 < a < 2 in the following 
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analysis, and the integration in Eqs. m and m yields the fractional Laplace 
operator for the FSE Q in the form of the Riemann-Liouville fractional 
derivative [17 1 

cL 






2T(2 - a) cos ^ 


di 


\x-y\ °‘i^iy)dy. 


( 12 ) 


Note that 


1 


2V{v) cos ^ 


(j){z)dz 
T - -r I 


is the Riesz fractional integral on the finite interval [a, b] I?! with a < x < b 
and 0 < < 1. It can be presented as the sum of the left and right Riemann- 

Liouville fractional integrals 


(j){z)dz 


4>{z)dz 


Ja {x-zy- J, {z-xy-’'- 
2.1. Eigenvalue problem for the fractional Laplace operator 
Let us consider the eigenvalue problem 


= E-^ 


E 


(13) 


with boundary conditions e{x = ±L) = 0 that yields the solution of FSE 
([9]) tp{x,t) = J2E^Ee~^^*'i’E, where coefficients qe are defined by the initial 
conditions ipoix) = 'if’{x,t = 0). We rewrite the fractional Laplace operator in 
the form FSE ([91), which is convenient in the following analysis, 

i-L r 


/ E 1 rOO 


f{y)dy. (14) 


Here we use again that on the real axis |A:p = kf. One easily finds that the an¬ 
tisymmetric (odd) eigenfunction can be found in the form (see Appendix) 




TO = 1,2 . 


(15) 


/L L 

which satisfies the boundary condition = ±L) = 0. Substituting solu¬ 

tion CSl) with f{y) = '^E{y) in Eq- (|14p and performing integration straight¬ 
forwardly, one obtains 

/ *- \ a 

miTX 


ra ■ ^EX 

L sm —-— = Dc 


/ hm-K \' 

) 


sin ■ 
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Therefore, is the eigenfunction with corresponding eigenvalue 


E°‘^^ = = D, 


odd 


/ hrmr \ ' 


- 


(16) 


The same procedure is performed to find symmetric (even) eigenfunctions 


»T,even _ »T,even 

— ^2m+l 


: COS ■ 


(2m + l)7r 


with the eigenvalue 




Vl 

h(2m + l)7r 


2L 


2L 


(17) 


(18) 


f. fl, 


2.1.1. A comment on the ground state 

This result coincides with the solutions obtained in Ref. [Ij, but the crucial 
difference here is that solutions (USD - (0 belong to a completely different 
fractional operator (ED, which is defined on the finite-size range of the potential 
well. Note also that in this case a deficiency with the ground state cos 


correctly stated in Refs. [2|, no longer exists, since it is obtained by the 
appropriately determined operator. 


3. Fractional path integral: fractional quantum mechanics with topo¬ 
logical constraints 

One should bear in mind that the Feynman-Kac formula in Eq. (|TD and FSE 
([SD are introduced in unbounded Euclidian (configuration) space. In contrast, 
the configuration space of a particle in the infinite potential well is contracted 
to the finite size of the potential well. This leads to the quantum mechanics in 
a space with topological constraints. Even for the conventional local quantum 
mechanics, the path integral presentation is not an easy task, as stated in [l^ . 
Fortunately, since eigenvalue solutions (ED-(ED are known, the path integral 
for the Levy process in the finite area can be constructed. Therefore, one can 
obtain Eq. © with infinite limits of integration for the Riesz operator. Let us 
define for convenience |a;) = ^'(((^‘^(a:). Then, Eq. o can be rewritten: 

V’(x, At) = {x\ip{M)) = = J dxi{x\e~"^"^*^^\xi){xi\i;o) . 

(19) 
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We focus on the evolution Green’s function 


G{x,At;xi) = {x\e-^^^^*/^\xi)= [ {x\y). 

J-oo J-L 

( 20 ) 

Taking into account that 

^ CSO 

(x\y) = ^ ^ sin(fc„x) sin(fc„j/) = — (21) 

n—0 

and the Poisson summation formula 

OO OO 

E E S{x-m), (22) 


1 — — C 


m=—OO 


one obtains (see also Q) 
Gix,At;x,) = - E / 


— ih^" ^\k\^DaAt 


^ik{x—xi-\-2Ll) _ ^ik{x-\-xi-\-2Ll) 


1 ^ 

lE E 


dk 


— ih"^ 


Z — -^X 1 — — 00 ‘ 

X exp {ik{z — xi + 2LI) + Z7r[6*(—z) — 0(a;i)]} 


(23) 


where the 0 functions in the exponential provide the correct signs for z = ±x. 
Now, we take into account 


1 °° 

jE E 




(24) 


which returns one to the fractional integration with infinite limits, as in Eq. O- 
The essential difference here is the appearance of the topological phase due the 
9 functions in Green’s function (1231) . Continuing this procedure by repeating it 
N times {t = NAt), one arrives at the correct analogy of Eq. ([3]) for the infinite 
well potential 


G{x,t;xo) = T,Z-c 


.E.=±.+2^lUUI- 


dxi 

OO J 


n iv+i coo 
J-c 


dkj 


X exp{f [kjixj - Xj-i) + TT9{-Xj) 


■n9{—Xj-i) — h°‘ ^\kj\°‘DaAt^^\‘^b) 
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where xat+i = z. The topological term is a purely boundary expressioio, and 
therefore, integration over all Xj in Eq. (1^ yields a product of 6 functions 
S{kj—kj+i), which finally yields the Fourier inversion of the Green function 
in the form of Fox’s function UQ in its Fourier form 

" 1 


G{x,t;xo) = ^ 


1 — — C 


47r 


POO 

^ik{x—xo-\-2lL) ^—ik{x-\-xo-\-2lL) 

/ —oo 

■ 


(26) 

Because of the properties of Fox’s function, at a = 2, the Green function 
reduces to the free particle in the box with the Green function [l^ 


1 — — 00 *- 


. {x-XQ+2lLy -im (x + XQ+2lLy 

t — 6 * 


(27) 


4. Example: Adiabatic approximation in moving walls 

Handling the path integral expression in such a simple form as Eq. (IMll . 
one can consider a system of a Levy particle inside moving walls. For 0 = 2, 
this problem corresponds to the so-called Fermi acceleration [l^ , where chaotic 
dynamics can be realized because of the interaction of nonlinear resonances. 
However, for 1 < o < 2, the physical implementation of the classical limit of 
the problem is vague, since a particle is not free. Therefore, the interaction 
of nonlinear resonances is not considered. Here, the problem is treated in the 
adiabatic approximation. In this case, Eq. (12611 can be used for further quantum 
mechanical analysis. Let the infinite walls at a; = ±L move periodically, such 
that the boundary conditions for the wave function are il){x = ±L(t)) = 0, 
where L{t) = L + e with e/L <C 1. To calculate the density of states, we 

^For a = 2, this expression coincides with the results in Ref. |lfil| . For 

the self-contained presentation, we shall adjust some formulae thereof. Namely, 
taking the continuous time limit, one arrives at the path integral G{x,t; xq) = 
h E“-oo T.z=±:r+ 21 L I Dx(t) f exp {|S}, where the action 5[d = drfpx - 

/f(|p|) — h7rxS(x)] contains the new topological term ^topf-s] = —hn dT[x5{x)] = h7r[0{—z) — 
0(—fco)], which is a purely boundary expression. 













are interested in the trace of the Green function 


g{t) = J dxoG{xo,t-,xo) 

where G{x = Xg, t; a:o) can be obtained from Eq. (I^Sl) in the form 


G{xo,t;xo)= ^ 


1 — — C 


47r 


nOO 

^ik{2lL) ^—ik(2xo-\-2lL) 

y —OO 

■ 


Performing summation over I by the Poisson summation formula 


E 


^^2iklL 


- i + ^E' 


2iklL 


11: 


(28) 


(29) 


(30) 


1——00 1—1 m——oo 

one can perform integration over k, which leads to summation over the spectrum 
k = km = TTm/L. Then, performing integration over xg, one obtains for the 


trace 

OO 

ga{t) = 

m—1 

One also obtains this result from the expression for the Green function G{xg^ t; xg) 

exp(—ii?°‘^'^t/S)4'°‘^‘^(a;o)4'*°‘^‘^(xo). The Fourier transform go(^') = ■^['?o(i)] 

yields the density of states pg{E) = —^^g{E) 


pg{E) = hYHE-Dc,\hkmn 


(32) 


where E = E°‘^^, which satisfies the eigenvalue problem, discussed in Sec. 12.11 
In the case of the moving walls, at the replacement L —^ L{t), the adiabatic 
approximation makes it possible to obtain the trace of the Green’s function in 
the form 

g,{t) = = Y , (33) 

m 

where we take h= 1. Let us present a shift operator in the form 


^(Esin vt)- 


= E 


n= —OO 


h 




(34) 


Then, the Fourier transform over time can be easily performed, which yields the 
following density of states 


Pe{E) = 


OO 


27T 


^ 2LYHE - nv - Da 


TTITI 


‘) • (35) 


n= —OO 
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Parameter^ determines the quasi-energy spectrum, namely, its band. Therefore, 
for the moving walls the discrete spectrum (fT51) changes essentially: 


Em ^ Em,n{^) — Tll^-\- Dq 


7Tm 


L -h £ sin ^ 


: nv-\- Do 


nm 


Trm 

— saDa ^ , sm ^ . 




(36) 


This solution corresponds to the reconstruction of the energy spectrum of the 
stationary problem m to the quasienergy spectrum with a narrow band struc¬ 
ture. 


5. Conclusion 


The introduction of the quantum Levy process in a box leads to the need to 
account for the topological constraints in the space. This problem can be treated 
in the framework of the path integral formalism with the Levy measure. A 
correct path integral consideration is possible when the eigenvalue problem p3p 
is appropriately defined for the infinite well potential (|S]), and the eigenfunctions 
are known, or can be obtained for the fractional (nonlocal) operator (fl^ . An 
analytical expression for the evolution operator is obtained in the path integral 
presentation. For a = 2, the path integral has the limit of the correct expression, 
which is well known in local quantum mechanics with topological constraints 
[l^ . Although the results of the eigenvalue problem coincide with the solutions 
obtained in Ref. [li , the solutions (fTSl) - ( H]) belong to the fractional operator 
which differ from those obtained in [l[ , since these are defined on the finite- 
size range of the potential well and correspond to the Levy walks. In this case, 
the ground state cos is correctly defined as well, and to some extent, this 
ends the discussion on the ground states (see Refs. 0, [sj). 

An important point of any “fractional quantum mechanics” is the classical 
limit fi —>■ 0, where a serious question emerges about the physical meaning of its 
classical counterpart. Obviously, a quantum Levy particle is not a free particle in 
the classical limit. As an example of such a possible realization, the Levy process 
is considered for oscillating walls. For a = 2, this problem corresponds to the 
so-called Fermi acceleration [l^ leading to classical and quantum chaos. But 


10 











for 1 < a < 2, the physical implementation of the classical limit of the problem 
is vague, since a particle is not free. Therefore, the interaction of nonlinear 
resonances is not considered here. The problem is treated in the adiabatic 
perturbative approach. This question can be an interesting task for future 
studies of classical and quantum chaos of the Fermi acceleration mechanism in 
the nonlocal quantum mechanics. 

This research was supported by the Israel Science Foundation (ISF). 


Appendix A. Eigenvalue problem 


Let us consider the eigenvalue problem (fT51) for the fractional Laplace oper¬ 


ator and find the antisymmetric (odd) eigenfunctions 'Lgin the form 




m = 1,2 ,... 


(A. 1) 


VL L 

which satisfies the boundary condition = ±L) = 0. We rewrite the 

fractional Laplace operator in the form FSE ([9]), which is convenient in the 
analysis, 

rL r 




>-L 






'fTiy)dy. 


Let us first perform integration over y. This yields 

,2'Km sin(A:L) 

LB - ■ -e “ j = 

~L 


1 

■‘-I r ikyA-iz iky—izl / i \r 

^ ] = (-!) 


(A. 2) 


(A. 3) 


iL {k + z){k — z) 

where z = irmy/L. The next step is integration over k. From Eqs. (lA. 21) and 
(|A. 31) . we have integrals 
1 



ik(L-x) 

-i(L+x) 

/_oo {k + z){k-z) 




= /(+) - /(-) 


(A. 4) 


where sign (-I-) corresponds to the analytical continuation in the upper half 
plain, while (—) corresponds to the analytical continuation in the lower half 
plain. Using the Residue theorem, one obtains that the integration yields 

(A. 5) 


1 _ 1 prtQ 1 

fmiT \ 

\ L ) ' 

k17^} 
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Acting on this result by the rest part of the operator, which reads (^) {idx), 

one obtains 







(A. 6) 
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